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ABSTRACT 

Electromagnetic field together with zero-mass charges moving in this field form a well-behaved 
semi-dissipative dynamical system - Electrodynamics of Massless Charges (EMC). We give equa- 
tions of EMC, argue that EMC is an adequate theory for calculating pulsar magnetospheres, give 
an illustrative numerical calculation (showing that bolometric luminosity of an aligned rotator 
is approximately equal to half the spin-down power) . EMC looks like a portion of the full pul- 
sar theory that will resolve the already calculated bolometric luminosity into light curves and 
spectra. 



1. Electrodynamics of Massless Charges 
(EMC) 

EMC postulates that positive and negative 
charges move with unit, meaning speed of light, 
velocities: 

E x B ± {B Q B + E E) 

v± = WTeI ' (1) 

Here Eq and Bq are the proper electric field scalar 
and the proper magnetic field pseudo-scalar, de- 
fined by 

Bl-El = B 2 -E 2 , B E Q = B • E, E o >0. (2) 

Equations ([2| mean that E and |i?o| are the field 
magnitudes in the frame where the fields are par- 
allel, with the sign of Bq negative when the fields 
are antiparallel. 

The " Aristotelean" dynamics ([!]) , with the par- 
ticle velocity rather than acceleration given by the 
force, results from an infinite radiation damping of 
massless charges. Equation ([I]) means that in the 
frame where the fields are parallel, charges move 
along the common direction of the fields at the 
speed of light. Moving along the fields nullifies 
the leading-order radiation damping term. And, 
because of curvature radiation, the parallel dy- 
namics is also instantaneous, fixing the velocities 
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The full EMC system of equations is, in the 3+1 
split, 

B = -V x E, (3) 

E = V x B - p + v + + p_v_, (4) 

P± + V • (p±v±) = Q. (5) 

Here p± are the particle densities (multiplied by 
the absolute value of the charge) , and the source Q 
represents pair creation/annihilation. The source 
Q can be chosen at will. This is because our mass- 
less particles do not carry energy and momentum 
(again because of the strong radiation damping), 
and therefore pairs can be created and annihilated 
arbitrarily. The only requirement is that Q > if 
one of the densities p± vanishes - there can be no 
negative particle densities. 

Somewhat surprisingly, even in this incomplete 
form, EMC seems to be useful, because for a broad 
range of assumed sources Q, the EMC predictions 
remain virtually unchanged, and one can indeed 
choose Q as he wishes. In pulsar applications, an 



lr The terminal Lorentz factor due to curvature radiation is 
given by 7 4 ~ EoR 2 /e, where R is the radius of curvature 
of the trajectory. Our massless charges are actually mov- 
ing slower than light. But since 7 is given by the ratio of 
macroscopic and microscopic quantities, in all possible ap- 
plications of EMC one has 7 > 1, allowing to approximate 
the actual velocity by Jxl. 



arbitrarily chosen Q allows to calculate only the 
bolometric luminosity, while for the spectra and 
light curves one must calculate the source Q from 
the actual charge-photon kinetics. This means 
that the photon distribution function needs to be 
added to the EMC system. The important point, 
however, is that even in this full pulsar theory, 
the field and charge dynamics is still given by the 
EMC system p|5 |. 

Applicability of EMC to pulsars requires in- 
stantaneous parallel dynamics at light cylinder, 
which gives 



L > L e 



R 



2/3 



(6) 




where L is the spin-down luminosity, R is the light 
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cylinder radius, r e — is the classical electron 
radius, L e = ^ = 8.7 x 10 16 erg/s is the "classi- 
cal electron luminosity" . Equation (JsJ) , for a 30ms 
pulsar, gives I> 6 x 10 30 erg/s, and is easily sat- 
isfied. 

2. Aligned rotator in EMC 

To illustrate how EMC works in pulsar applica- 
tions (the only application the author is aware of) , 
we simply assume that Q, whatever it is, creates 
a prefixed quantity of high-multiplicity plasma, 
meaning 

P++ P- = \p\+ fPo, (7) 

Where p = p + — p_ is the charge (Goldreich- 
Julian) density, / is the multiplicity and po > \p\ is 
a fiducial charge density greater than the absolute 
value of the actual charge density \p\. 

In the figure we show the results for / = 10 
and po = \JB 2 + E 2 /R : where R is the spheri- 
cal radius. This po simply represents the charge 
density which is capable of changing the fields by 
an order unity factor. The simulation technique 
is similar to what we have used when simulating 
pulsars by Strong-Field Electrodynamics (SFE) 
(Gruzinov 2011a). The light cylinder is at twice 
the stellar radius. 

We get an absolute replica of the pulsar magne- 
tosphere calculated in the high- conductivity limit 
of SFE. This calls for an explanation. 



Fig. 1. — Isolines of (i) the poloidal magnetic flux 
- thick black, (ii) the electric potential - thinner 
red, (iii) the poloidal current - thin blue (shown 
only below the equator for clarity of the picture). 



3. EMC and SFE 

SFE is Maxwell theory plus Lorentz-covariant 
Ohm's law (Gruzinov 2011a and references therein) 



J = 



pE x B + 



j 2 a 2 El(B B + E B) 
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(9) 



Here a is the conductivity, which in SFE is as- 
sumed to be an arbitrary scalar depending on Eq 
and Bq. This Ohm's law has been postulated just 
from the Lorentz covariance requirement. 

SFE can be derived from EMC if one postulates 
the Lorentz-invariant source term 



Q = k{g 2 E% - 2#j„- 



(10) 



where j± are the Lorentz-covariant 4-currents 
j± = p±(l,v±). In the limit k — > oo, in the 
frame where the fields are parallel and the charge 
density vanishes, the particle densities become 
p± = ctEq/2 and we do get the current 

It is good to know that SFE (a simpler but 
postulated theory) is derivable from EMC (a first- 
principle theory, but with more fields and not 
fully defined). This does not explain, however, 
why an arbitrarily chosen high-multiplicity plasma 
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of §2 gives the same magnetosphere as the high- 
conductivity SFE. 

The universality of the high-conductivity SFE 
magnetosphere follows from the expression for the 
electric current in EMC: 



pE x B + ( P+ + g-KgpB + £ E) 
5 2 + El 



(11) 



Now comparing the EMC current ( 11 ) to the SFE 
current we see that the equilibrium state of 
the high-multiplicity plasma of §2 is identical to 
the SFE equilibrium with a certain (weird) space- 
dependent conductivity - but as long as the con- 
ductivity is everywhere high, its actual value is 
irrelevant. 

The following seems to be a fair analogy. The 
temperature jump at the ideal gas shock can be 
calculated in dissipative gas dynamics with arbi- 
trary but small viscosity. Similarly, the amount of 
damping of an ideal pulsar magnetosphere can be 
calculated in EMC with arbitrary but high multi- 
plicity, or in SFE with arbitrary but high conduc- 
tivity. 

EMC and SFE give identical results for the 
bolometric luminosity. But only EMC seems to be 
a portion of the full pulsar theory, because EMC 
gives a faithful description of the particle density 
(except for creation and annihilation if any) . 



layer is weaker in the inclined case (Spitkovsky 
2006). 

I thank Matt Kleban for useful discussions. 
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4. Summary 

EMC - Electrodynamics of Massless Charges - 
is a well-behaved dynamical system. The charges 
carry zero energy-momentum and obey Aristote- 
lean dynamics. The system is semi-dissipative. 
The electromagnetic energy can decrease, but for 
certain field configurations it is precisely con- 
served. 

For pulsars, in the high-multiplicity limit, EMC 
gives an undamped Poynting flux everywhere ex- 
cept at the equatorial current layer. As the 
charges carry no energy, the equatorial Poynting 
flux damping must be interpreted as bolometric 
luminosity. For an aligned rotator, one gets the 
bolometric luminosity ~ 50% of the spin-down 
power (Gruzinov 2011b). For a generic inclina- 
tion angle, the bolometric luminosity has not yet 
been calculated. Presumably the bolometric lumi- 
nosity at generic inclination is smaller than 50% of 
the spin-down power, because the singular current 
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